Suppression of Spontaneous Supercurrents in a Chiral p-Wave Superconductor 
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The superconducting state of SRO is widely believed to have chiral p-wave order that breaks 
time reversal symmetry. Such a state is expected to have a spontaneous magnetization, both at 
sample edges and at domain walls between regions of different chirality. Indeed, muon spin resonance 
experiments are interpreted as evidence of spontaneous magnetization due to domain walls or defects 
in the bulk. However, recent magnetic microscopy experiments place upper limits on the magnetic 
fields at the sample edge and surface which are as much as two orders of magnitude smaller than the 
fields predicted theoretically for a somewhat idealized chiral p-wave superconductor. We investigate 
the effects on the spontaneous supercurrents and magnetization of rough and pair breaking surfaces 
for a range of parameters within a Ginzburg-Landau formalism. The effects of competing orders 
nucleated at the surface are also considered. We find the conditions under which the edge currents 
are significantly reduced while leaving the bulk domain wall currents intact, are quite limited. The 
implications for interpreting the existing body of experimental results on superconducting SRO 
within a chiral p-wave model are discussed. 
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I. INTRODUCTION 

Strontium ruthenate, Sr2Ru04 (SRO), has attracted 
considerable experimental and theoretical study since 
its discoveryJi It was the first perovskite supercon- 
ductor to be discovered which did not contain cop- 
per and is believed to have unconventional pairing 
symmetry. ■^'^ Numerous experimental results have been 
interpreted as evidence of a superconducting order pa- 
rameter with spin-triplet pairing^i^i^i^ii and broken time 
reversal sy mmet ryi^i^ i ^^^^^ The simplest order parameter 
consistent with these observations corresponds to a chiral 
p-wave, Px ± ipy, Cooper pairing symmetry;^ analogous 
to the A phase of superfluid He-3.^^ 

This state is expected to give rise to spontaneous 
supercurrents flowing along the sample edgC j^^i^'^'^^i^^ 
which are screened by the Meissner effect so that the 
magnetic field is zero inside the superconductor. The 
net result is a magnetic field confined near the edge of 
the sample. These spontaneous currents and fields can 
also occur within the sample at domain walls between 
Px + ipy and px — ipy domains. ^S-iH Surface currents and 
domain wall currents in a chiral p-wave superconduc- 
tor have been studied by Matsumoto and Sigrist^!^ and 
other a"'^^'"'^^i"'^'''i"'^^i"'^^'^° and should be observable by scan- 
ning probe measurements 1^ As well, muon spin reso- 
nance experiments have been interpreted as evidence for 
internal fields present at domain wallsi^i^. 

Recent scanning Hall bar and superconducting quan- 
tum interference device (SQUID) microscopy measure- 
ments did not see the expected signatures of spontaneous 
currents at the sample edges and surfaces <22i2^ These null 
measurements set upper limits on the spontaneous cur- 
rents which are approximately two orders of magnitude 
smaller than the values predicted from simple chiral p- 
wave order 3^ Given the considerable body of experimen- 
tal results taken as evidence for chiral p-wave order, it is 
important to understand whether the absence of observ- 
able magnetization at the edges can be explained within 



a theory of bulk chiral p-wave superconductivity. One 
possibility discussed by Kirtley et al^ is domains at the 
surface smaller than 1 or 2 microns on average. Given 
the size of the experimental probes, this could account for 
the null resultsi^ Indeed, Josephson tunnehng measure- 
ments were interpreted as evidence of chiral p-wave order 
with small dynamic domains, although other results 
would be incompatible with such small domains at the 
surfaceiiii or in the bulkii The formation of domain walls 
is energetically unfavorable in the Meissner stat o^^'^° and 
the samples are considered clean (otherwise Tc is notice- 
ably reduced as expected for unconventional pairing^!), 
so such small domains arising from dynamics and pin- 
ning would be somewhat surprising. However, an alter- 
native to the Meissner state, one which favors domains 
of roughly the size of the penetration depth, has been 
proposedi^ Additional experiments are required to ei- 
ther rule out or confirm the presence of small domains. 

Alternatively, one might expect surface roughness or 
other surface effects to reduce the spontaneous currents 
and, in this paper, we investigate this possibility. Pre- 
viously, only ideal (specular) surfaces of a chiral p-wave 
superconductor have been considered although the 
effect of a rough surface has been considered for a neutral 
chiral p-wave superfluid^^ where screening currents are 
absent. Rough surfaces can be studied in the Bogoliubov- 
de Gennes (BdG) formalism or closely related Greens 
function formalism. Here, we use a Ginzburg-Landau 
(GL) formalism allowing us to more readily study the ef- 
fect of a variety of surfaces as the parameters in the the- 
ory are varied. These correspond to studying different 
microscopic Hamiltonians in the BdG formalism, which 
each stabilize a Px i ipy superconductor. The BdG for- 
malism is more accurate at low temperatures, although 
for specular surfaces it was found that the GL calcu- 
lations gave qualitatively similar results for the sponta- 
neous currents and fields. ^^'^^ We also consider the effect 
of surfaces which nucleate a non-chiral p-wave order pa- 
rameter, while maintaining px ± ipy in the bulk, as a 
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possible mechanism for suppressing the predicted edge mensionless variables takes the formic 
currents. 



II. GINZBURG-LANDAU EQUATIONS 

The Ginzburg-Landau free energy functional describ- 
ing a single layer of Sr2Ru04 expressed in terms of di- 
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in 
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Here, u and v are the x and y-components of the order 
parameter, respectively. The bi and ki are dimensionless 
material dependent constants, and di = — iai are 
the usual gauge covariant derivatives. The position x is 
scaled by the coherence length, f , a is the dimensionless 
vector potential, and we have introduced k — the 
usual GL parameter. Parameters satisfying &2 > and 
63 — 4^2 < stabilize the chiral p-wave state. 

The values of the coefficients in the free energy can be 
computed in the weak coupling limit of a BCS supercon- 
ductor with triplet pairing aligned along z as described 
by Furusaki et ali^ and correspond to 62 = |, 63 = 0, 

— \: k2 — \t kj, = J and ^4 = i. It is also sometimes 
convenient to introduce the variables k± — 5(^3 ± k^). 



We take A = 190 nm and ^ = 66 nm as parameters appro- 
priate for strontium ruthenate, unless noted otherwise. 

The order parameters are parametrized by u — \u\e'^^ 
and V — \v\e^'^^^'^\ We require the free energy to be sta- 
tionary with respect to variations of the order parameters 
and the vector potential to obtain 6 coupled non-linear 
partial differential equations. We consider the case of a 
boundary at x = with the superconductor occupying 
the half-plane x > 0. The symmetry in the problem al- 
lows us to discard y derivatives, as well as to choose the 
gauge where Qx — 0. Therefore, u and v are taken to 
be functions of x only. This reduces the problem to the 
solution of the following 5 equations: 



= - fci|u|" - {k+ + {ay\v\ sin(0))' - + (i + 62 - \h) \u\' + + 62 cos(2<^) + \h) \v\M 

+ ki\u\9 + k2\u\a'l^ ~ k+ cos{4))\v\ay {26 + ^ — fc+ sin((/))ay |w| -f fc_ cos((/))|w|ay0 + /c_ sin((/))ay |w| , (2) 

= - fc2b|" + {k+ - fc-) {ay\u\ sin(0))' - + (1 + 62 - i&3) \v? + + 62 cos(20) + \h) \u\M 

+ ki\v\a'^ + k2\v\{9 + 4>) — k+cos{(f>)\u\ay {26 ^ + /c+ sin((/))ay|u| + k- cos{(j))\u\ay(l) -I- fc_ sin(0)aj^|u| 



(3) 



(4) 



= - fc2 {\v\^{e + cj))'^ + (fc+ - fc_) {\u\\v\ayCos{(j))) - 2|up|i;p62sin(0)cos((/)) + fc+|u||w|ay sin(0) {26 + 
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I 

We integrate Eq. [5] to obtain: which is then used in Eqs. [^Hl] to reduce the problem 

/,M II I I I 12 1' to only 4 equations. These equations are solved self- 

, ^ 2k+ cos{4>)\u\\v\ay - k2\vY(j) ^ ^ 

fcilwp + felwp ' ^ ' 
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consistently using a numerical relaxation algorithm sim- 
ilar to that described by Thuneberg,2I The equation for 
the current follows from Eqn. [5] and the Maxwell equa- 
tion VxB = — j. We identify the terms in the current 
proportional to the vector potential as the screening cur- 
rents, and the others as the spontaneous surface currents. 



III. WEAK COUPLING RESULTS 

We consider a superconductor filing the half plane 
X > with a surface at a; = 0. To derive the bound- 
ary conditions on the order parameters, we follow Am- 
begaokar et alM By considering quasiparticle trajecto- 
ries, they show that a surface is always pair-breaking for 
the component of the order parameter which is normal 
to the surfacei^ This implies the boundary condition, 
m(0) — 0. Combining this with the restriction that no 
current should pass through the interface leads to the 
condition 



\v\ 



= const. 



(8) 



x=0 



For a clean surface (specular scattering) next to an insu- 
lator the appropriate choice is const. = Oj^ 

The boundary condition on |m| suppresses it near the 
surface, and the term (| — ^2 + jb^j \u\'^\v\ in the equa- 
tion of motion for \v\ allows \v\ to vary close to the sur- 
face. In general, if — a) = [j — b2 + jbs) > 0, \v\ wiU 
be enhanced near the surface, and if — a) < 0, it will 
be suppressed For weak coupling, this coeSicient takes 
the value | and, as seen in Fig. (TJ the y component of 
the order parameter, is larger at the surface than in 
the bulk. 

It is also noteworthy that even if we start from an ar- 
bitrary relative phase, (p, between the two order param- 
eter components, self consistent solution of the bound- 
ary problem forces the relative phase to ±^ even as the 
magnitude of the order parameters vary spatially at the 
surface. That is, the Px ± ipy state is maintained near 
the boundary. 

An analysis of the GL equations for fixed </> = f shows 
that there are only two conditions under which the spon- 
taneous current vanishesj^ (i) 4^2 — b^ = Q and (ii) 
ki = ^2, fcs = ^4, which gives |w| — \u\ everywhere. 
The first case is a boundary of stability on the chiral p- 
wave state which requires 4^2 — &3 > 0. Thus, it follows 
that only with |m| — \v\ everywhere can the current be 
reduced to zero while maintaining chiral p-wave order in 
the bulk. Since by symmetry |m| is suppressed at the sur- 
face, we must introduce an effect which also suppresses 
\v\ to reduce the current. As a candidate, we examine 
the effects of a rough surface with variations on a scale 
much smaller than the coherence length, so that it can be 
treated as a boundary condition. For such a surface the 
treatment of Ambegaokar et. al still applies^^, u(0) = 0, 
and \v\ satisfies Eqn. [51 
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FIG. 1: Self consistent solution of the GL equations for the 
weak coupling parameters, (a) The x and y components of the 
order parameters scaled by the bulk order parameter. Here 
the subscripts s,d and pb denote the case of specular, diffuse, 
and pair breaking scattering, respectively, (b) The magnetic 
field and current distributions scaled by and -s^^tt re- 
spectively. A comparison of the integral of the magnetic field 
over 25 coherence lengths shows a 22% reduction for the dif- 
fuse case compared to specular scattering. 



Following deGennes^^ the constant in Eqn. [8] is de- 
noted as -J , where b — oo corresponds to specular scatter- 
ing and the limit of diffuse scattering provides a minimum 
value of 6 = 0.54^^ We also consider an extra suppres- 
sion of the order parameter corresponding to 5 < 0.54 
which could be caused by magnetic scattering at the sur- 
face which disrupts the triplet pairing. The limiting case 
6 = corresponds to a completely pair breaking surface 
with both components of the order parameter driven to 
zero at the surface. 

The result of a self-consistent solution of the GL equa- 
tions for the weak coupling parameters is shown in Fig. 
[T] for both specular and diffuse scattering as well as 



for the fully pair breaking boundary condition. The x- 
component of the order parameter is almost the same 
in all three cases since any surface along y is fully pair 
breaking for this component. As x approaches the sur- 
face, the y-component of the order parameter still grows 
up as the x-component is suppressed, but it also is ul- 
timately suppressed close to the surface due to the pair 
breaking boundary condition. This behaviour can be at- 
tributed to the different healing lengths of the two com- 
ponents in response to a perturbation in a;. These quali- 
tative shapes of the order parameters replicate those from 
previous work on the effect of a clean surface^ and of 
a rough surface on a neutral chiral p-wave superfluid^SS, 
This demonstrates that the GL theory and boundary con- 
ditions treated here are in good agreement with the mi- 
croscopic BdG and Green's function calculations. 
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IV. RESULTS AWAY FROM WEAK COUPLING 

Since the limit 4^2 — &3 — > causes the currents to 
vanish, we examine self-consistently the dependence of 
the solutions on the parameter 452 — 63- In Fig. [2] we 
demonstrate the dependance of the fields and current for 
a value of the parameter 4^2 — 63 which is closer to zero 
than for the weak coupling parameters. We first notice 
that the healing length of the parameters is extended in 
this regime, resulting in a broader spontaneous current 
distribution near the edge. As this parameter is tuned 
closer to zero the spontaneous fields become progressively 
smaller until the state is no longer stable and they vanish. 
However, the rough surface boundary condition has a 
smaller effect on the reduction of the fields. A reduction 
in magnetic field by changing parameters in this way will 
reduce all the magnetic signatures, and will not be able 
to account for the experiments taken as evidence for time 
reversal symmetry at domain walls. 

Variation of the coefficients fc^, which are the stiffnesses 
of the order parameters, also has an effect on the mag- 
netic fields produced at the surface . The weak cou- 
pling ratio 1^ = 3 is apparent in both Figs. [l]and[2l as 
the x and y components heal over different length scales. 
In Fig. [3] we change the parameter k\ and observe the 
change in integrated magnetic field as the two order pa- 
rameters are forced to change on the same length scale. 
Allowing fci and k-2, to become equal does not change 
the integrated magnetic field significantly if the bound- 
ary conditions on u and v differ, as they do for specular 
or diffuse scattering. On the other hand, if the surface 
is pair breaking and v is significantly suppressed at the 
surface (in addition to u) then the integrated field falls 
off much faster. For fci ~ k2 the currents and fields are 
zero. 

It is interesting to ask if the parameter choice ki — k2 
will also have a large effect on the currents produced 
at a domain wall. There are two types of domain 
walls whereiiS, (I) the relative phase continuously changes 
through the wall or (II) the phase changes discontinu- 




FIG. 2: Self consistent currents, fields and corresponding or- 
der parameters for the parameters &2 = 63 = |, fci = |, 
k+ = J. In this parameter regime the currents and fields are 
naturally suppressed, with the integrated magnetic field 23% 
less than that of weak coupling parameters. Also, the changes 
in the currents and fields due to surface roughness is reduced 
here resulting in only 10% change in the integrated magnetic 
field. 



ously. Recently a paper studied domain walls in the GL 
formalism and found a parameterization which connected 
the two types of domain walls 3* They expressed the mag- 
netization as a function of the GL parameters. Previ- 
ously, Matsumoto and Sigrist showed that the domain 
wall configuration of type II is energetically favoured.''^ 
Since only one order parameter is driven to zero for this 
type of domain wall, our previous analysis shows that a 
non-zero magnetization will be produced, consistent with 
the result of Logoboy and Sonini^ This situation would 
allow for the magnetic signals attributed to domains in 
the bulk, as well as for the lack of currents at the edge. 
The difference in magnitude of ki and fc2 is associated 
with the different energy costs of longitudinal and trans- 
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FIG. 3: Dependence of the integrated magnetic field on the 
parameter ki — k2 for specular and diffuse scattering, and the 
pair breaking boundary condition, all other parameters are as 
for weak coupling. The weak coupling parameters for fci, k2 
correspond to fci — A;2 = | . 



FIG. 4: Computed currents and fields for a solution with spa- 
tially varying relative phase near the surface. The parameters 
used are those for weak coupling, with 62 negative for 2; < 5 
but positive elsewhere. 



verse perturbations respectively. There there is no sym- 
metry which would require that fci = fc2. In general one 
would expect the longitudinal fluctuations to be stifFer, 
as is the case for the weak coupling parameters. 

One last feature of Fig. [3] which stands out, is the 
increase of the overall magnetic signal for the specular 
boundary condition as fci ^ fc2. To understand this we 
express the screening currents (js) in terms of the sum 
and difference of fci and k2'. 

k'^Js = -ay{ki\v\'^ + fc2|Mp) 

= -^[(H^ + i-n + (fci-fc2)(H^-i«n], (9) 

where we have used the constraint ki -\- k2 = 1. From 
Eq. [5] we see that as fci — fc2 — > the screening cur- 
rents are reduced in the region where (jvp — jitp) > 0, 
which is only satisfied near the sample edge. The change 
in fci — fc2 also changes the spatial dependence of the 
spontaneous currents. The spontaneous current both in- 
creases in magnitude and is pulled closer to the sample 
edge as fci is reduced. This move of the spontaneous 
currents to the region where the screening currents are 
reduced results in an increased magnetic signal. 

V. EFFECT OF COMPETING SURFACE 
ORDER 

To consider the effect on the edge currents of compet- 
ing order which is favoured by the surface, we consider 
several possibilities. First, we allow the parameters that 
stabilize the chiral p-wave state to vary spatially near the 
surface. In particular we change the sign of the 62 term 
in the free energy in a region near the surface (62 > 



favours the Px + ipy state, and 62 < the Px + Py state). 
The Px +Py ground state does not support spontaneously 
generated supercurrents, and being stabilized near the 
edge could reduce magnetic signatures. The resulting 
currents and fields from a self-consistent calculation are 
shown in Fig. [4l Notice that both the overall magnitude 
of the fields is suppressed, and there is a change in sign. 
This alternating magnetic field is much like at a domain 
wall and will also affect the magnitude of the measured 
magnetic fields as described by Kirtley et alM^. While 
these spatially alternating magnetic fields could reduce 
the measured signals, this scenario has a different order 
at the surface and is incompatible with the interpretation 
of some of the surface tunneling measurements 

We also consider the effect of subdominant non-chiral 
p-wave order coexisting with the chiral p-wave state. 
This can be modeled by adding the following terms to 
the free energy: 

/2 = a2\w\^ + b^\w\^ + fcs {\Dxw\^ + \Dyw\''^ 

+ b5\w\^ + \v\^) + be [w*^ {u^ + V^) + C.C.] . 

(10) 

These terms would be caused by the addition of any one 
of the other unitary states allowed under the crystal sym- 
metry. However, symmetry does not allow gradient terms 
coupling the new order parameters to the old ones. This 
means that even when the parameters are such that this 
new order parameter can grow up near the surface, it has 
little effect on the shape of the old order parameters and 
hence, causes insignificant changes to the currents and 
fields, unless the transition temperature for this compet- 
ing order (determined by a2)is high enough. 

Since we wish to maintain chiral p-wave order in the 
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bulk, we consider an 02 which varies spatially. In particu- 
lar we consider the case where the of w is greater than 
that of the chiral p-wave near the surface. Self consistent 
solutions of the GL equations show that the new order 
parameter grows up at the surface, suppressing the chiral 
p-wave state which is recovered in the bulk. This con- 
figuration gives rise to an alternating magnetic field, of 
similar magnitude to that shown in Fig. [51 with the max- 
imum magnitude of the fields being about 30% smaller. 

VI. DISCUSSION 

In an attempt to reconcile the results of various ex- 
periments on strontium ruthenate, we have explored a 
variety of mechanisms that could be responsible for re- 
ducing the spontaneous edge supercurrents generated by 
a + iPy superconductor, while simultaneously main- 
taining the px + ipy state in the bulk. One can think of 
several possibilities for reducing all spontaneous currents, 
at edges and at domain walls, such as multiband effects 
or using GL parameters near the boundary of stability for 
Px + 'i-Py However, one would then need to look for al- 
ternative explanations for the /xSR measurements which 
have been taken as evidence for chiral p-wave domain 
walls. Therefore we have focused on effects which reduce 
the edge currents, but not the domain wall currents. 

In particular, we examined the effect on the sponta- 
neous supercurrents of rough and pair breaking surfaces, 
as well as the dependence on the GL parameters. The 
spontaneous edge currents are zero or vanishingly small 
only in a very small region of parameter space in the 
presence of a fully pair breaking surface. This tuning of 
parameters is unlikely to be realized in a physical system 
as it requires the coefficients for longitudinal and trans- 
verse gradients to be equal. 

The effect of nucleating a non-chiral order parameter 



at the surface, while maintaining chiral p-wave order in 
the bulk, was also investigated. This can give rise to so- 
lutions where the magnetic field alternates in sign near 
the surface. These alternating magnetic signatures could 
produce null results for the edge currents if the length 
scale of the alternating magnetic field was sufficiently 
short. Again, this scenario would be difficult to reconcile 
with the tunneling measurements.^'^'' 

Leggett has proposed an alternative wavefunction that 
reduces to the BCS wavefunction for the case of s-wave 
pairing^!.. While, for the chiral p-wave case, the BCS 
wavefunction predicts an angular momentum of the con- 
densate of Cooper pairs given by i^^i^i^ Leggett finds 
that the angular momentum of the condensate described 

by his wavefunction is ^ (^) ' "^^^^ large suppression 
of the angular momentum would indeed reduce sponta- 
neous surface currents. However, this would also result in 
a suppression of all magnetic signatures and thus would 
leave the positive /iSR results unexplained. 

In summary, we have identified a number of possible 
explanations for the absence of observable edge currents 
in a chiral p-wave superconductor. However, each of 
these is then inconsistent with the interpretation of tun- 
neling and/or /iSR results that are interpreted as a direct 
observation of the fields induced by supercurrents at do- 
main walls. It would be illuminating to investigate this in 
greater detail by modeling the /iSR lineshapes expected 
for chiral p-wave domain walls. 
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